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Decaying topological defects, in particular cosmic strings, can produce a significant flux of high 
energy neutrinos, photons and cosmic rays. According to the prevailing understanding of cosmic 
string dynamics in an expanding Universe, the network of long strings loses its energy first into 
string loops, which in turn give off most of their energy as gravitational radiation. Ifowever, it has 
been suggested by Vincent et al. (VHS) that particle emission may be the dominant energy loss 
channel for the long string network. In this case, the predicted flux of high energy particles would 
be much larger. Here we calculate the predicted flux of high energy gamma rays, neutrinos and 
I cosmic ray antiprotons and protons as a function of the scale of symmetry breaking rj at which the 

. strings are produced and as a function of the initial energy mj of the particle jets which result from 

' the string decay. Assuming the validity of the VHS scenario, we find that due to the interactions 

with the cosmic radiation backgrounds all fluxes but the neutrino flux are suppressed at the highest 
energies. This indicates that the observed events above the GZK cutoff can only be accounted for 
in this scenario if the primary particle is a neutrino and rj is somewhat less than the GUT scale, 
i.e. rj < 10^^ eV. The domain of parameter space corresponding to GUT-scale symmetry breaking 
, is excluded also by the current observations below the GZK cutoff. A new aspect of this work is the 

CO . calculation of the spectrum of the tau neutrinos directly produced in the decay of the X-particles. 

This significantly increases the tau neutrino signal at high energies in all cosmic string scenarios. 

CN \ 

K*- ! PACS numbers: 98.80.Cq, 98.70.Sa, 98.70.Rz, 98.70.Vc, 96.40.Tv 

O . I. INTRODUCTION 

00 

' Measurements of the spectra of high energy 7— rays, neutrinos and cosmic rays have emerged as a useful constraint 
. on particle physics theories which predict topological defects (for recent summaries see e.g. Refs. 0, II). Although 
II ■ it does not appear likely that the emission of high energy particles by topological defects can explain the observed 
Oh spectra, the current data can be used as upper bounds to constrain theories of particle physics beyond the Standard 
^ Model. 

Of particular interest for cosmology are theories giving rise to cosmic strings. Some time ago, the spectrum of cosmic 
rays from non-superconducting strings was computed |4[ |, |^ under the assumption of the standard scaling picture 
10, ^, H] emerging from studies of cosmic string dynamics in an expanding Universe (see Refs. (l0[ |ll|, |l^ for recent 
reviews). In the standard picture, the long string network evolves into string loops which then decay predominantly 
' by gravitational radiation. It was found that for strings with Gji ^ 10"^ (required for strings to be relevant for cosmic 
structure formation p3| , p^), where /i is the mass per unit length in the string and G is the Newton gravitational 
constant, the predicted particle fluxes are substantially lower than the observational detections or limits. 

Recently, however, Vincent, Hindmarsh and Sakellariadou have challenged the standard picture of string evo- 
lution. They claim that the small-scale structure on the strings does not scale with the expansion of the Universe. 
In the VHS scenario, the long strings lose their energy directly into particles instead of string loops. This leads to a 
greater production rate of particles and hence to the expectation of greater fluxes of 7— rays, neutrinos, and cosmic 
rays at Earth. 

Expanding on these points further, detailed numerical simulations |l9| have demonstrated that in an ex- 

panding Universe, the network of long cosmic strings {long meaning with curvature radius greater than the Hubble 
radius) approaches a scaling solution in which the number of long string segments crossing each Hubble radius ap- 
proaches a constant value v. In the VHS scenario, the string network maintains this scaling solution by the direct 
radiation of scalar and gauge particles. In the standard scenario, the string network maintains the scaling solution 
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by continuously giving off part of its energy in the form of string loops with radius smaller than the Hubble radius. 
Based on certain theoretical arguments (see e.g. it is expected that the distribution of loops in the standard 

scenario also eventually takes on a scaling solution, i.e. when all lengths are scaled to the Hubble radius, the statis- 
tical properties of the distribution of string loops are independent of time. However, the resolution of string network 
simulations (Nambu action strings) is not yet good enough to be able to verify the standard scenario. Moreover, 
recent field theory simulations ||2l[| provide some evidence that the string defects give off their energy not in the form 
of string loops, but directly by scalar and gauge particle radiation, thus corroborating the VHS scenario (These latter 
simulations remain controversial [p2| , p3|). 

In the standard cosmic string scenario, the string loops oscillate due to the relativistic tension, and decay slowly 
by the emission of gravitational radiation. Only a small fraction of the energy is released in the standard scenario in 
the form of scalar and gauge particle radiation. This occurs via the process of cusp annihilation , during the final 
loop collapse or by the evaporation of black holes created from string loops |^ . It has been shown ||, |[ ^ that 
for GUT-scale (G/i ~ 10~^) strings the flux of ultra-high energy (UHE) particles in the standard picture from these 
mechanisms is at or below the relevant observational measurements or limits. The predicted flux increases, and may 
be observable, if the symmetry breaking scale is significantly lower. 

In the VHS scenario, in contrast, all of the string energy is released directly from the long strings as scalar and 
gauge particle radiation. Additionally, the predicted flux increases as the symmetry breaking scale increases. Hence, 
the flux of high energy 7— rays, neutrinos and cosmic rays in the VHS scenario is expected to be much larger than 
the fluxes from loops in the standard scenario. Indeed, it has been remarked ^ that the fluxes from GUT 

scale strings in the VHS scenario will be above the observational limits. In this paper, we point out that this limit 
comes only from the neutrino emission at high energy and the cascades produced by the electromagnetic emission, 
whilst the nearest expected VHS string is too distant to produce an observable flux of cosmic rays or 7-rays. In the 
VHS scenario the highest energy fluxes (above the GZK cutoff) of cosmic rays and 7-rays are suppressed due to the 
interaction of these particles with the photons of the cosmic radiation backgrounds. The result is that no isotropic 
flux of cosmic rays and 7-rays from a network of VHS strings is expected to be observed. However, the neutrino 
emission from low G/i VHS strings may produce an observable signal. We present a detailed calculation in this paper 
of the neutrino, photon and cosmic ray fluxes for a range of string and decay scales. We also show that the predictions 
depend sensitively on the initial energy mj characterizing the decay of the particle radiation, which in turn depends 
on the presently unknown physics in the regime between the scale of electroweak symmetry breaking and the scale of 
string formation. We find that decreasing mj leads to a decrease in the upper cutoff in the predicted spectrum and 
an increase in the cosmic ray fluxes at intermediate energies. Thus, even the additional freedom of decreasing mj 
cannot make GUT strings obeying VHS dynamics consistent with the observational constraints. Nevertheless, in the 
event that a segment of long string is substantially closer than the expected average distance to the nearest string, a 
highly anisotropic flux of cosmic rays protons, antiprotons and 7-rays could be observed above the GZK cutoff in the 
VHS scenario. This possibility, however, is not supported by the current observations of events above the GZK cutoff 
which exhibit large scale isotropy and at most only small scale clustering in arrival angle |28| . 

In this paper we perform for the first time the calculation of the tau neutrinos directly produced by the hadronic 
decays of the X-particles. The calculation of the direct Vt spectrum has been not previously been included in the 
analysis of any cosmic string scenario. This was because it was assumed that the Vt spectrum was highly suppressed at 
all energies compared with the and spectra. Instead only the much weaker v-r spectrum generated by the cascade 
of the X-decay i^e and off the cosmic relic neutrino background was described. Recently, however, MacGibbon, 
Wichoski, and Webber |3^ have shown that, in hadronic jets at accelerator energies, the tau neutrinos are a 
significant fraction of the total neutrino spectrum at high x, where x is the ratio of the neutrino energy to the 
jet energy. Extrapolating their fragmentation function parametrization for Vr to UHE energies, we show here that, 
in the VHS and all other cosmic string scenarios, the directly produced tau neutrinos generate a signal at Earth 
comparable at the highest energies with the expected Ve and signal. This flux is orders of magnitude greater than 
the contribution from the cascade Vr ■ One consequence of our result is that the detection of a significant fraction of 
Vt in UHE neutrino events may be due to hadronic decay at the source and not —^ Vr oscillation in transit. 

In the following sections, we present the neutrino, photon and cosmic ray spectra for vanilla emission taking — 
and omitting SUSY particles in the jet decay. However, we discuss in Section IV the effects of including a SUSY 
sector, non-zero neutrino mass, neutrino oscillation or other extensions to the Standard Model. In the Appendix, we 
perform the full computation of the flux of secondary muon and electron neutrinos resulting from muon decay in the 
particle jets. We use units in which ?i c = G = 1, unless otherwise stated, and take a Hubble constant of ho = 0.65. 
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II. JET FORMATION AND FRAGMENTATION 

If the network of long strings scales, then the corresponding energy density poo at time t is 

Poo = vpt^^ , (1) 

where fi is the mass per unit length in the string. This equation defines the constant v which counts the number of 
strings per volume in the scaling solution. The value of v can be determined from simulations of cosmic string 
network evolution. Current work |l^ gives v ~ 13. In the VHS scenario, the scaling solution is maintained 

not by the production of string loops, but by the radiation of gauge and scalar particles (which we collectively call 
X-particles in the following). Making use of the equation of state of a string, the continuity equation becomes 

where nx and mx are the number density and mass, respectively, of the X-particles. mx is effectively the symmetry 
breaking scale at which the string formed, 

mjf ~ ry - ^1/2 (3) 

(e.g. mx 10"'^^ GeV for Gp ^ 10~^). Strictly speaking, p^^'^/rj ^ 0{1) with the exact factor depending on the form 
of the potential of the scalar field. 

The decay of the high energy X-particles will lead to the production of jets analogous to the QCD jets seen in 
accelerators. To simplify the discussion we assume that all jets in the X-decay have the same initial energy, mj. In 
this case, the decay of a single X-particle will lead to mx/mj jets, and the number density of jets generated by the 
energy release of long strings is 

at 6 mj 

The largest uncertainty in the calculation of the flux of high energy particles generated by strings comes from our 
ignorance of the structure of the jets at ultra-high initial energies. Jet production has been studied in detail in QCD. 
From such studies, the fragmentation function of jets into photons, neutrinos and baryons is known |3^ , at least 
to a good approximation, up to a few TeV. However, in our case the particles producing the jets are superheavy 
scalar and gauge particles. Following Refs. |3^, the fragmentation functions of such jets are usually calculated 
by extrapolating the QCD fragmentation functions to higher energies (in which case mj is taken to be mx/2). This 
is appropriate if there is no new physics which does not match this extrapolation between the electroweak symmetry 
breaking scale and the unification scale (the scale at which the defects are produced), but it is not justified if there is 
new physics which does not sufficiently match the extrapolation. We will parametrize ignorance about new physics 
and its effects on the jets from the X-particle decay by introducing a new scale nij (which may be the scale of the new 
physics) as the scale to which we can satisfactorily extrapolate the fragmentation functions from the QCD regime. 
We will then assume that the initial scalar or gauge X-particle produces mx /mj jets with initial energy mj each. 
(In reality, there may be intermediate decay steps between the mx and mj decay scales.) While this may or may not 
be the precise description of the decay, it will at least help further quantify the dependence of the final fluxes on the 
uncertainties in the fragmentation process. 

The initial jet particle decays into quarks and leptons on a time scale of fimj^, where a is the coupling constant 
associated with the physics at the energy scale mj. The quarks then hadronize on a strong interaction time scale. 
More than 90% of the total energy in a QCD jet goes into pions, with the majority of the remainder going into 
baryons (mainly neutrons, antineutrons, protons and antiprotons). On astrophysical time scales, the neutral pions 
decay into two photons, the charged pions decay into electron and muon neutrinos and electrons or positrons, and the 
neutrons decay into protons and leptons. In addition, a relatively small but still significant number of tau neutrinos 
are produced by the decay of heavy quarks and tau leptons. Non-jet leptons Ix produced by the decay X — > q + lx 
may also influence the observable 7-ray flux below E ^ 10^^ eV (see Section |l|). 

The distribution of energies E of the primary QCD jet decay products (predominantly equal numbers of 7r'^,7r^ 
and TT^) can be approximated by the fragmentation function |34] 

where x = E/mj is the fraction of the jet energy carried by the decay product and E continues down to ~ 10^ 
eV. The two body decays of the primary tt", tt"*" or tt" decay particles lead, on the average, to two photons and two 
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primary muon neutrinos for every three pions. The energy distribution of the decay products can be obtained by 
integrating ^ from ax to 1 with the invariant measure dx/x, 



dN 
dx 



5 .16 

— rr 



2a 



1/2^1/2 



ia 



'1/2^-1/2 



(6) 



where the constant a depends on the decay process being considered (see below). Eq. (||) apphes to the final spectrum 
of photons produced in the jet. The photons result from tt" decay. In this case, both of the decay particles have 
vanishing rest mass, and it is possible for a single photon to carry away the entire pion energy in the lab frame. Hence, 
a = I for the photons. 

Electron and muon neutrinos result from the tt"'" and 7r~ decay. In the decay process tt^ fj,^ 
products (the muon or antimuon) has nonvanishing rest mass. In this case, the integration limits when integrating 
are not x and 1 (see e.g. and The correction to the upper integration end for i/g and i/^ is negligible 

but not so the change in the lower integration limit for i'^ which becomes ax with 



I'll, one of the decay 



a 



2.34. 



(7) 



when i?7r ^ ~ 140 MeV, the charged pion mass, and ~ 106 MeV is the muon mass. Equation ^ describes 
the and neutrinos immediately produced by the decay of the charged pions. The final spectrum of v^, and i/^ 
neutrinos, however, is comprised of the i/^ and Pp, neutrinos produced by the subsequent decay of the fi~^ and /i^, 
together with the contribution from those immediately produced by the decay of the 7r+ and 7r~ . The final spectrum 
of Ve and 9^ neutrinos is produced by the decay of the muons and antimuons. Taking into consideration the spin 
polarization of the muons, we show in the Appendix that the final spectra of the Ve and neutrinos generated by 
jet decay are approximately 
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for the muon neutrinos produced in the first stage of the pion decay 
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for the muon neutrinos produced by the decay of muons in the pion decay, and 
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for the final electron neutrinos. 

Equation (||) also applies to the primary baryonic decay products such as protons and neutrons. Assuming that 
3% of the total energy of the jet goes into baryons and antibaryons and noting that most of the energy of decaying 
neutrons and antineutrons is transferred to daughter protons and antiprotons, the final distribution of protons and 
antiprotons is approximately 



dx 16 



(11) 
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FIG. 1: Comparison of the three jet fragmentation functions discussed in the text. The solid line represents the fragmentation 
function (|^) used in this paper to calculate all ultra-high energy particle fluxes but v-r , the dash-dotted line is the from of @, 
the dotted line is the approximation (^3|) of (^2|), and the dashed line is (p^. 



Equation (|^) is one approximation to the numerically computed QCD fragmentation functions, taken from Ref. p4[ . 
It is derived simply by requiring that the number of decay products scales as E^^^ and applying energy conservation 
and scale invariance. Further the ratio of nucleon to pions is assumed at all x to equal the ratio of total nucleons to 
total pions per jet. Even with these reservations, we find that Eqs. (|^) - match the i/,, and spectra generated 



in 10 TeV e ' e 
10-^ <x <1 
Another form is 




q q events simulated by the HERWIG Monte Carlo code ||3^ to within a factor of 2 in the range 
There are other analytical approximations of the fragmentation functions used in the literature. 



ax 



which can be well approximated by 
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3.322x 
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(12) 



(13) 



for values of x between 10 < x < 10 ^. A third formula is based on the Modified Leading Logarithmic Approxi- 
mation (MLLA) [pSl of QCD. At smaU x, the expression becomes [|9l Eol 



where 



dN' ^ K r In'^x/xn 
dx X ^ 2a^ 



„ 2 1, /"^J^3/2 



(14) 



6"'^' A 

Xm — (A/mj)^/^ and A — 0.234 GeV. Conservation of total energy in the jet 

1 dN' , 
X— — dx = 1 

dx 

gives the normalization constant K. In all approximations, extrapolations and Monte Carlo simulations of the 
fragmentation functions at high energies, the fraction of the jet energy carried by nucleons (or equivalently the 
nucleon multiplicity) lies in the range 3-10% but the precise value and dependence on x and mj is unknown p4|. 
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Fig. 1^ shows a comparison of the three fragmentation functions (|^), (|T^ and (^4|). The difference between 
these functions is small in the energy range 10~^ < a; < 10~°'^, and wiU be unimportant for order of magnitude 
considerations. Hence, in the following we will use mostly approximation (H). We stress however that no extrapolated 
fragmentation function should be regarded as accurate to more than an order of magnitude at highest and lowest x. 
Investment in precision in those regions of the predicted spectra is unjustified and misleading. 

For the first time in cosmic string emission analysis, we also include the distribution for the tau neutrinos produced 
in the jet decays, which had erroneously been assumed to be negligible in all previous work. While the number of 
produced per mj > 1 TeV jet is less than 10~^ of the number of and f^, the high energy Vr are predominantly 
produced by the initial decays of the heavier quarks with shorter lifetimes and the h'e and v^, are produced by the final 
state cluster decays of the much lighter pions. This leads to significantly greater relative contribution from Vt at high 
X than previously assumed. The fragmentation distribution (|^) is no longer relevant for the tau neutrino. Instead 
MacGibbon, Wichoski, and Webber [E9, 30| find that the fragmentation function for production in 300 GeV - 100 



TeV e'^e ^ qq events, as simulated by HERWIG |32[| can be parametrized as 
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Consistent results are obtained from jet decay events simulated by PYTHIA/JETSET ||3l| . 

The expressions for the number density of jets (P, and for the energy distribution of the jet decay products, (|^), 
®; (|fO|); ( pr]) or (15), can be convolved to obtain the expected flux F{E) of high energy photons, neutrinos, and 
cosmic rays with energy E from a VHS string distribution 



F{E) = 



in 



tc 



, , dn J , , ,. 
dt'^m 



-3 dN dE'[E,z(t')] 
dE' dE 



(16) 



In Eq. (|l^), particles observed today {t = ip) with energy E were produced at a time t' with energy E' = E'[E, z{t')], 
where z{t') denotes the redshift at time t' . The factor (z(t') + l)^^ expresses the dilution of the particle number density 
in an expanding Universe. The lower cutoff time tc{E) corresponds to the maximal redshift from which particles of 
present-day energy E can reach us. This cutoff can be due either to interactions with the ambient cxtragalactic and 
galactic media during propagation, or to the constraint that the initial energy E' must be less than the initial jet 
energy mj. The upper cutoff time tistr corresponds to the latest time that the emission from the VHS cosmic string 
network can be considered isotropic (see Section [II). Thus our upper cutoff time is that from which emission will 
have just isotropized by today. 

For protons, in order to calculate the Jacobian '^'^ ^dE^*"^^ apply the continuous energy loss approximation (GEL) 
and solve Eq. (^6|) numerically. The use of the GEL approximation is acceptable in this case where the distance to the 
sources is much larger than the attenuation length for the particle in the cosmic radiation backgrounds (see Section 

For photons, if electromagnetic cascades are not taken yet into account (see below), we obtain (see Eq. (g) with 
a = 1) 
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For neutrinos, if cascade off the relic neutrino background is not included, we have (see Eq, 

15 f E \ ^ 



47r z \iJij / 



+ 1)) 

E \-3/2 



0.24 {{Zmax + 1)"'/' - {Zrmn + 1)"'/') f ) 

2.98 ((Z„,a. + 1)'/' - {Zrmn + 1) ) f — ) ^'^^ 

1.38 ((z™,, + 1)3/2 _ (^^^^^ ^ 1)3/2) I' AV^' 

\ nij / 

— I 

mj / 



+ 0.14 {{Z„,ax + If - {Z,mn + if) ( j 

- 0.03 ((Z„,a. + if - {Z,mn + lf)( — Y 

\mj/ 



E \2 



for electron neutrinos and (see Eqs. (H) and M 
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for muon neutrinos. 
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Similarly for the tau neutrinos, the flux from the string distribution neglecting cascading is (see Eq. (h 
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In Eqs. (17-201), Zmin is the redshift corresponding to tistr a,nd Zj^axi^) is the redshift corresponding to tc 



TTl J 

Zrnax{E) + 1 = r77m(Zco(£^) + 1, ^) , (21) 

where Zco{E) denotes the redshift cutoff due to astrophysical interactions. For convenience, we have presented Eqs. 
(0 - 1^) for the c&se Zfjidx 

[E) < Zeq. In our figures, we have calculated the more general result for Zjnax{E) > Zeq 

where appropriate. 

The first important conclusion to draw from Eqs. (p"7]-p^) is that for energies E substantially smaller than the 
effective initial jet energy 777,,/, the fragmentation term proportional to {E/mj)~^^'^ similarly dominates the flux of 
photons, electron and muon neutrinos giving 

E'F^,,^+,^^,^+,^{E) (X 7777^/' (22) 

In the case of the tau neutrinos the fragmentation term proportional to (E/mj)^^^'^ dominates the flux so that Eq. 
(|^) gives, 

E^F.^+c,^ cx 7777'/' . (23) 

Hence, for flxed cosmic string mass per unit length ^, the spectrum of high energy photons and neutrinos increases 

as the jet energy r77j decreases. This occurs because for smaller values of 777,,/ there are more jets emitted from the 

3/2 

string and the fragmentation function as a function of E is more compressed. These two effects overcome the mj 
1/2 

or 777.y factor of the dominant term of the fragmentation function. We also note that the fluxes in the VHS scenario 
increase as G/i increases, or equivalently the symmetry breaking scale increases, in contrast to the general decrease 
seen in the standard cosmic string scenario. 

To complete the calculation of the photon, neutrino and cosmic rays we need to know the values of Zmax{E). For 
protons and antiprotons with energy between '--^ 5 x 10^^ eV and ~ 5 x 10^^ eV, the dominant interaction is e+e" pair 
production off the cosmic microwave background. At higher energies, photopion production becomes the dominant 
energy loss mechanism. The onset of photopion production corresponds to the GZK cutoff energy |Q . Above this 
energy, the rate of energy loss increases dramatically, and nucleons do not reach us from cosmological distances. 

Applying the results of |Q, we find, in the case of protons and antiprotons, Zcop(E) for various values of 777 j = 77/2 
(where -q is the scale of the symmetry-breaking) corresponding to G/i = 10^®, 10^*, lO^^*', 10^^^^ 10^^^, and 10^^® as 
shown in Fig. ||. 

In Fig. ^ we plot the average maximum distance that a ultra- high energy proton or antiproton can travel in the 
CMB. After propagating ^ 100 Mpc the proton or antiproton energy decreases below ^ 6 x 10"'^^ eV irrespective of 
the initial jet energy. For an expanded treatment see e.g. 

For photons of energies greater than about 10^^ eV, Zco^{E) is determined by pair production off the cosmic photon 
background. The resulting function Zco^{E) is shown in Fig. |^ (taken from Ref. [^). For lower energies, pair 
production off nuclei determines the form of Zco ■ 

For neutrinos, the dominant process is scattering off the relic 1.9°K cosmic background neutrinos. In this case, 
Zco„ {E) is given by Q 

/ E\-^'^ 

Zco„ =i 9 X lO^f — j , E<'iy.lQ^^eV (24) 
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FIG. 2: The average maximal redshift 



, (E) from which protons and antiprotons can reach the Earth as a function of the 



energy at arrival. The solid lines represent from top to bottom G/i =10 ,10 
decrease of the curve above 4 x 10^® eV corresponding to the GZK cutoff. 
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FIG. 3: Energy of a proton (antiproton) as a function of the distance of travel through the CMB, for the same 6 values of Gfi 
(determining the initial energy) used in the previous figure. 
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(25) 



The decay of lepton pairs produced in the neutrino scattering process further enhances the spectrum of all neutrino 
species and allows neutrino emission, corresponding to present day energies E < 10^^ eV, to be detected from 
significantly higher redshifts than those given in (|2j) and (^5|) [Q. The cascade enhancement due to decay is only 
relevant at high redshift where the neutrino mean free path is sufficiently small. In the case of electron and muon 
neutrinos, we will neglect the cascade enhancement of Zco^ because the dominant term in the neutrino flux from VHS 
strings, (|l^) and (p^, is proportional to (1 — {zmax,, + and the modiflcations to the v^, u^, and Vr spectra 
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FIG. 4: The maximal redshift Zco^{E) from which photons can reach the Earth as a function of the energy at arrival 

expected at Earth are negligible above E ^ 10^'^ mx [Q- The cascade enhancement at low energies increases as mx 
increases due to the greater energy in the initial jets. 

III. COMPARISON WITH OBSERVATIONS 



As we mentioned above, the diffuse fluxes of particles in the VHS scenario are produced directly by the decay of 
X-particles. The X-particles decay immediately after they are radiated from the long string segments. The result is 
that the particle fluxes are produced along the strings. The number of long strings per Hubble volume v we noted in 
Section I is approximately 13. As the Universe expands the comoving intcr-string distance grows as 



Dints 



:(1 



-1/2 



(26) 



Fig. 1^ depicts the evolution of the average distance between long string segments. The distance between the 
observer and a long string segment, Dis is estimated to be 



D,, 



If the distance today between the observer and the long string segment is larger than the maximum average propagation 
distance for particles arriving with an energy E then the fluxes at Earth are exponentially suppressed |^ . Based on 
Eq. (pfl) we see that presently the average distance to a long string segment is 



0,ls 



2000 Mpc, 



(27) 



which corresponds to a redshift zig ~ 0.58. 

The present average redshift to a long string segment zis can be taken as a rough estimate of Zmim i-e., of the 
redshift corresponding to the latest time (tistr) that the emission from the VHS string network can be considered 
isotropic. This sets an upper limit of the integral (|l^) and consequently in Eqs. (pT|-po|) we have 



Zls 



0.58 . 



(28) 



Hence only the particle fluxes emitted from the VHS network before Ustr [z > Zmin) are relevant for the diffuse flux. 

A segment of long string could be substantially closer than Do, is, the expected average distance to the nearest 
string. In this case, a flux of cosmic rays protons (antiprotons) and photons above the GZK cutoff might be observed. 
However, these fluxes would be highly anisotropic in contradiction to the present observations of UHE events above 
the GZK cutoff p8i. 
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FIG. 5: Average distance between long string segments as a function of the redshift. 

A. Protons 

In what follows, we discuss the proton flux. Because the fragmentation functions and relevant interactions with 
the cosmic backgrounds are to first order charge symmetric, the expected antiproton flux from the VHS strings is the 
same as the proton flux. 

The average maximum distance at the present time to the source of a proton arriving with energy ~ 4 x 10^^ eV 
(the onset of the GZK cutoff) is < 160 Mpc (see Fig. ||). This is equivalent to the source being at a redshift z < 0.036, 
much smaller than the distance from which the flux can be considered isotropic (~ 2000 Mpc). It implies that no 
isotropic flux of protons from VHS strings is expected at Earth above the GZK cutoff. On the other hand, a diffuse 
proton flux from the strings should reach the Earth at arrival energies i? ^ 6 x 10^^ eV, well below the UHE event 
energies. Although no isotropic flux of protons at UHE energies is expected at Earth, the UHE protons emitted at high 
redshifts will have interacted with the universal radio background (URB), CMB, and the infra-red/optical (IR/0) 
backgrounds to produce secondary fluxes of 7-rays, electrons, and neutrinos. These secondary fluxes could in principle 
be used to constrain the VHS model. However, in the context of the TD scenarios including the VHS scenario, these 
nucleon-induced secondary fluxes are negligible compared to the corresponding fluxes directly produced by the decay 
of the X-particles. This is because in the X-particle decay, far more direct 7-rays and neutrinos are produced than 
nucleons. 



B. Photons and Charged Leptons 

The maximum redshift at which a photon arriving today with an energy E could have been emitted is given in 
Fig. ^ If we compare this redshift to the present average distance to a long string (z/s ~ 0.58) we would deduce 
that the isotropic UHE 7-ray flux should also be exponentially suppressed. In the case of 7-rays, this is not the full 
story. Fig. |^ describes the interaction length of the photons, i.e., the distance traveled by a UHE photon before it 
is absorbed in the cosmic backgrounds. The UHE photons initiate electromagnetic cascades as they travel in the 
cosmic backgrounds. An important consequence is that the effective penetration length of the cascade is considerably 
larger than the interaction length In the VHS scenario, both the interaction length and the effective penetration 
length are much smaller than the distance to the nearest expected string (cf. Eq. (^)). The result is that the 
E > 10^^ eV 7-ray are essentially subtracted from the arriving flux: The flux above E ^ 10^° eV is depleted due 
to interactions with the URB; the flux in the range 10"^^ ^ E < 10^° eV is depleted by the CMB; and the flux in 
the range 10" < E < lO^'^ eV is depleted by the IR/0 backgrounds. The energy lost by the photons emitted with 
energies above ~ 10^^ eV at high redshifts is recycled into energies below E > 10^^ eV. Analytic calculations |Q have 
shown that the electromagnetic cascade spectrum has a generic shape below 10^^ eV, i.e., the secondary cascade 
radiation spectrum is insensitive to the injection spectrum. This implies that, even in the absence of the detection 
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of direct UHE photons, the 7-ray flux below E > 10^^ eV can be used to constrain the VHS scenario. The charged 
leptons created in the X-particle decay also initiate electromagnetic cascades on the cosmic backgrounds. As in the 
case of the UHE photons, the energy of the UHE charged leptons is recycled and enhances the photon flux below 
> 10" eV. 

The constraints on the electromagnetic emission in TD models with a spectral index p ~ 1, which is applicable to 
the VHS scenario, stem from the limits on the present diffuse 7-ray background between 3 x 10^ eV and 10^^ eV 
These require that |l|, ||, || 

Qem ^ 3 X 10^23 ^^-3 gg^-i ^ ^29) 

for h = 0.65 where Q%}^,i is the present rate of total energy injection into the electromagnetic channels. 

In the VHS scenario, we have from Eq. that the present rate of total energy injection into the electromagnetic 
channel (photon and charged lepton emission) is 

3 

eV cm"3 sec"^ . (30) 

Comparing to the Q%m limit in (^9[), we obtain a constraint on the energy density of a string in the VHS scenario of 

G^i < 3 X 10"" . (31) 

Here we have assumed an extragalactic magnetic field of B^xgai ^ 10~^^ G. Stronger values of the extragalactic 
magnetic field generate greater initial synchrotron losses at ultra-high energies, decreasing the UHE spectra, enhancing 
the cascade reprocessing into low energies, and thus lowering the upper limits on Q%j^,[ and Gp,. Because of the large 
distances to the sources, we have also neglected the uncertainties in the cosmic background URB and IR/0 intensities 
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C. Neutrinos 

The predicted flux of high energy neutrinos in the VHS scenario can be obtained from Eqs.(|l^ - |2^) using 
(p^), and ( |2^ ) to determine Zmax^- The resulting Ue and fluxes are shown in Fig. ^ and the flux in Fig. 
the case when mj = rj/2 where rj is the scale of symmetry breaking. 

The cascading of the tau neutrinos off the relic neutrino backgrounds will have little effect on the tau neutrino 
spectrum shown in Fig. |^. However, a secondary i/^- component peaking at lower energies will be created in collisions 
of the string-produced Vf, and i/^ with the relic neutrino backgrounds. This and i/^ cascading should increase 
the net flux below x ^ 10~^. The i^r signal expected at Earth from VHS strings will be the sum of these two 
components, as shown in Fig. ^ for Gfi ^ 10~^ |2^, |30j. In Fig. ^ we have plotted the cascade component calculated 
in Ref. | |4^ for a p = 1 TD scenario. This is only an estimate because the distance distribution of strings varies from 
model to model. We will address full modeling of the i^t fiux including cascades in a later paper. Note however that 
it is the high x region, where the cascades have little effect, which is the most observationally interesting. 

Our new result is that the Ve, i^fi and i>r fiuxes are of similar magnitude as £' — > to,/, unlike all previous string 
emission calculations which predicted a negligible Vr component at high x. In other string emission or X-decay 
scenarios, the Vf,, and v.,- flux should also be of comparable magnitude as a; ^ 1. Because of the uncertainty in 
extrapolating to high emission energies, the precise ratio of the neutrino fluxes at high x is not known (also, as plotted 
in Fig. ^ the extrapolation for Vr is derived from HERWIG Monte Carlo simulations while the and v^j^ functions 
are derived from the less rigorous approximation Eq. (^).) The actual ratio will be determined by the first step in 
the X-decay chain which is unknown and may be influenced by, for example, leptonic and SUSY decay channels (see 



Section IV) 



The most recent observational limits come from the EAS-TOP [g4|, the Frejus [q5|, and the Fly's Eye ||56| experi- 
ments. These experiments give upper limits on the cosmic ray neutrino flux in the energy range between 10^ eV and 
10^*^ eV. In Figs. ^ and ^ we have plotted the observational limits using the updated acci'f^N) cross-section By 
inspection, we see that the most stringent constraint comes from the highest energy flux in Fig. ^ which gives an 
upper bound on G/i of 

G^l < IQ-^". (32) 

This bound is derived at high x so it should be little affected by the cascades which are neglected in Fig. ^ 

We have also plotted on Figs. || and ^ the upper limits stemming from the planned sensitivities of the Pierre Auger 
project as estimated in Ref. pSf, and the OWL/Airwatch project as estimated in Ref. It is assumed that both 
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FIG. 6: Electron and muon neutrino fluxes in the VHS cosmic string scenario (for mj = r]/2) for various values of (lines 
from top to bottom, G^JL = 10"^, 10"*, 10"^°, 10"^^, 10"", 10"^^). The solid lines represent the v^, flux and the dashed lines 
represent the Ve flux. Points with arrows represent upper limits on the diffuse neutrino flux from the EAS-TOP [Q, the Frejus 
[ p5[ and the Fly's Eye |5^] experiments. The upper limits stemming from the projected sensitivities for both the Pierre Auger 
project and the OWL/Airwatch project ^] are also plotted. 




FIG. 7: Tan Neutrino flux in the VHS cosmic string scenario (for mj — rj/2) for various values of G/i (lines from top to bottom, 
G/i = 10"^ 10"*, 10"^°, 10"", 10"", 10"^*^). For each value of Gfi, the range of energies where the Ut from decays dominates 
the total i/t flux is depicted by the solid section of the line. Points with arrows represent upper limits on the diffuse neutrino 
flux as in the previous flgure. 



experiments will collect data over a period of a few years. Their sensitivities span the energy range between 10^^ eV 
and 10^^ eV and would allow more stringent upper limits on G/z up to 

Gfi < 10^12 ^ (33) 
in the case of null detection by the Pierre Auger project; and 

G/z < 10"" , (34) 
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in the OWL/Airwatch case. Unlike the present constraint (|3^), the Auger and OWL/Airwatch projects will constrain 
the I'f^ flux in Fig. ^ at energies E ~ 10^^ eV (i.e. the lower end of the energy range of these experiments). The 
Vr flux is below the present observational upper limits. The fr sensitivities of future experiments are expected to 
constrain G/i, but these constraints are expected to be an order of magnitude weaker than those from the i>e and 
fluxes. The slope of the i>r spectrum is steeper than the slope of the observational upper limits for both present and 
future data (cf. Fig. 0). On the other hand, the z^r component of a TD signature maybe easier to distinguish from 
other sources than the Ve and i/^ components. 

As is obvious from (|2^), ( p3| ) and Fig. ^the predicted flux of cosmic ray neutrinos increases as rnj decreases. Thus, 
the VHS scenario with GUT-scale strings cannot be made compatible with observational constraints by lowering the 
energy scale mj of the initial jets (which likely happens if new physics evolves between the electroweak symmetry 
breaking scale and the GUT scale). Raising mj above rj would alleviate the conflict, but this appears quite unnatural. 




FIG. 8: Neutrino flux in the VHS cosmic string scenario for various values of G/i — 10~ The solid line represent the jy^j flux, 
the dashed lines represent the Ve flux, and the dotted line the flux. The dash-dotted line represents the r^r from the cascade 
of Ve and due to the interactions with the relic neutrino background. Points with arrows represent upper limits on the 
diffuse neutrino flux as in the previous figures. 



The initial decay channels and branching ratios for the X-particle are unknown. The application of the extrapolated 
QCD-derived fragmentation functions to X-decay assumes implicitly that the branching ratios for X —^ qi + qi are the 
same as those for e~'"e~ — > Qi + qi. This may not be so for even hadronic decays at very high energies. The initial decay 
channel X q + Ix, where Ix is a charged lepton, should produce spectra which are decreased by no more than a 
factor of 2, compared with the X q + q spectra (assuming equi-partition of the X-particle energy between the initial 
q and Ix), with the exception of an Vr spectrum enhancement a.s E —> mx/'i, due to decay of the initial r ||2^ , 
The initial charged lepton Ix would also generate EM cascades which are effectively attenuated at high energies in 
the VHS scenario, as we remarked in Section II for the other EM decay products, but which would contribute to 
the spectra at low energies. This should increase the total energy going into the EM cascade channel by less than 
~ 50% (see Eq. (|30|)). Pure leptonic decay channels X ^ Ix + I'x^ ^ ^ + t^x or X ^ vx + i^'x would produce 
substantially suppressed spectra, because of the low multiplicity per decay compared with hadronic jets, except for 
a possible neutrino spectra increase as E ^ mx /2 due to the initial step in the decay chain. In the case of purely 
leptonic decays, the nucleon and 7-ray spectra would be wholly created by the collisions of the emitted leptons with 
the cosmic ambient media. It has been shown that in TD models generally, pure leptonic decays can only be made 
interesting if the cosmic relic neutrinos are massive, non-relativistic and locally clustered thus enhancing the relative 
cascade-generated component (Ref. [Q). However, it is highly unlikely that the X-particles (being high energy gauge 
and scalar particles) would decay purely into leptons in the initial step. Because of the much higher multiplicity of 
hadronic jets, any scenario containing reasonable initial branching ratios for the X-particle (say, to within a factor 
of 2-3: 30% pure hadronic, 30% pure leptonic and 30% hadronic and leptonic) would produce spectra dominated by 
the particles produced via the hadronic decay channels. Hence, even in the case of locally clustered relic neutrinos, 
the component produced in pure leptonic decays would be essentially dominated by the hadronic decays, with the 
possible exception of a UHE neutrino feature as E ^ mj5c/2. 
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FIG. 9: Electron and muon neutrino fluxes in the VHS cosmic string scenario for G/i = 10"^ and mj ~ lO^'^ , 10^^ , lO^S 10^'' , lO" 
eV (lines from bottom to top). The solid lines represent the z/^ flux and the dashed lines represent the flux. Points with 
arrows represent upper limits on the diffuse neutrino flux as in the previous figures. 



IV. EXTENSIONS TO THE MODEL 



In the analysis so far, we have assumed rrii, = 0. If neutrinos do possess a small mass, as present neutrino oscillation 
experiments hint, the predicted spectra from VHS strings are modified in two ways. Firstly, the interaction of the 
string- produced neutrinos with the 1.9°K relic background neutrinos is strongly enhanced at the resonance for 
production pOl , Eres — '^zl'^ "'^v ^ 4 x 10^^/(m,y/eV) eV. (If neutrinos are massless, the resonance occurs at 
-Eres ^ 3 X 10 eV [^.) The present upper limits from particle physics experiments on the neutrino masses are 
Wi^e ~ 3 eV, rriy^ < 0.19 MeV, and < 18.2 MeV The cosmological requirement that il„i < 1 places a 

stricter upper limit on the sum of the masses of the three neutrino species of < 92/i^ eV Taken together these 
limits imply that at least one resonance must occur for neutrinos today at some energy in the range 10^" eV 
^ Eres ^ 3 X 10^** eV. The exact value(s) of Ej-^s will depend on whether or not the neutrinos have mass, and the 
value(s) of the neutrino mass(es) if non-zero. If the neutrino masses are non-degenerate, the resonance will occur 
at two or three distinct i?res in this range. In the case that iJ^es < ^x, the resonance in the VHS scenario will 
produce a slight decrease in the predicted neutrino fluxes around Eres and in the energy decade below Eres, and 
an accompanying slight enhancement on all spectra below ~ \{i~^Eres due to secondary Z^ decay products. The 
amplitudes of all spectra should be modified by less than a factor of 2. (See for example Ref. |Q for the case of 
degenerate = 1 eV neutrinos in other cosmic string models.) If the relic neutrinos are non-relativistic and clustered 
locally on a scale less than the attenuation length for photons and nucleons, the effect on the photon, proton and 
antiproton spectra may be enhanced [6^ . 

Secondly, a non-zero neutrino mass may lead to oscillations between neutrino species as the flux propagates to 
Earth. In this case, the total number of neutrinos is preserved but the ratio of the flux in each species may depend 
on the neutrino energy, the distance traveled and neutrino masses. Recent observations by the SuperKamiokande, 
Kamiokande, MACRO and Soudan experiments of i/g and neutrinos produced by cosmic rays colliding with the 
Earth's atmosphere strongly suggest that neutrino oscillations (at least between the z/^ and Vr flavors) occur. [For a 
review of the current status of neutrino oscillation experiments and measurement implications see [ p4| and references 
therein]. The evidence for zenith-angle dependence in the SuperKamiokande data further strengthens this interpre- 
tation. Similarly, the long-standing observed deficit of solar neutrino flux may be resolved by i^e oscillation into v^^ 
or Vt although the solar neutrino experiments are not yet sufficiently consistent to favor a particular solution. The 
third set of evidence for neutrino oscillations comes from the LSND Collaboration accelerator experiment which has 
seen evidence for Vy, into oscillations. When combined with the atmospheric and solar data, this would require the 
existence of a fourth "sterile" neutrino species. The LSND measurements have yet to be independently confirmed. A 
number of new experiments are presently being commissioned with the objective of resolving these matters. The K2K, 
NuM-MINOS and CNGS experiments are expected to provide the definitive answer and exploration of the relevant 
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parameter space. In these experiments, neutrino beams will be sent from accelerators to detectors over baselines of 
hundreds of kilometers. K2K is already taking data and the NuM-MINOS and CNGS will be operational in 2002 and 
2005, respectively. 

Ignoring for now cosmological evolution, the probability of a relativistic neutrino undergoing vacuum flavor transi- 
tion when propagating a distance L is 



in the 3-flavor solution. The atmospheric data to date are best fit by an amplitude of sm^(6'o) = 0.82 — 1.0 and a 
mass difference of Am^^^ — (2 — 6) x lO^^eV^ |Q. For models with distant cosmological sources (e.g. TDs and 
AGNs) and a neutrino ratio at emission of F,^^ : F^^^ : F^^ ~ 1 : 2 : (in our case this occurs at << mj), the 
ratio of the spectra reaching the Earth will be F^^ : F^^ : F^^ = 1 : 1 : 1 in the 3-flavor model; similarly, in the 
4-flavor model consistent with the LSND data, an initial emission ratio of F^^ : F^^ : F^^ : F^^ ~ 1 : 2 : : 
will become F^^ : F^^ : Fy^ : Fy^ = 1 : 1 : 1 : at Earth |65). Since we showed earlier that the current average 
distance to the nearest long string is 2000 Mpc, even the highest energy neutrinos emitted in the present era will have 
undergone transition before reaching Earth. In this case, at << mj, both the and v-r flux reaching Earth should 
approximately equal half the i/^ flux expected without oscillations. In the case of the flux emitted in previous epochs, 
this will only be so provided the wavelength of the oscillation is much less than the mean free path of the and 
Vr in the Universe. If the oscillation wavelength is greater than the mean free path then the cascade structure will 
be determined by the original neutrino flavors at emission. The mean free path decreases sharply as z increases . 
For example, a lO^'^ eV neutrino at z ~ 10^ has a mean free path of ~ 100 kpc |Q. Since these too are the redshifts 
and energies at which cascading becomes relevant, cascade development should be predominantly determined by the 
neutrino flavors at emission. In the VHS scenario, however, the neutrino cascades have little influence on our regions 
of interest. Hence the predicted observable and v-r VHS fluxes in the case with neutrino oscillations should be 
approximately half the unoscillated flux at E << mj. As E ^ mj/2, i.e. in the region where the initial v,,, v^n 
and Vr VHS fluxes are of comparable magnitude, the ratio should also approach F^^ : F^^ : F^^ ~ 1 : 1 : 1 in the case 
with neutrino oscillations. The precise value of the modulated neutrino flux F^. as E ^ mj/2 will depend on the 
initial absolute and relative string-produced fluxes which we regard as uncertain by more than a factor of 2. 

Additionally, one possible solution to the solar neutrino deficit problem is vacuum v^. —^ or Vr oscillations with 
Auiq ^ lO^^^eV^ and sin'^26Q ^ 0.6 — 1.0. In this case, the phase of the transition probability, (j)Q = AtoqL/ {2hcE), 
is ^ 1 for i? > 10^^ eV neutrinos emitted today from a long string at 2000 Mpc. Thus the I'e spectra at Earth 
would be determined by Ve or Vr oscillation below E ~ 10^^ eV and exhibit negligible Vf. oscillation above 

E ^ eV. If this particular solar deficit solution is combined with the solution implied by the atmospheric data 
Am^(„ ~ 10"'^ eV^, an initial ratio at emission of F^^ : F^^ : F^^ ~ 1 : 2 : or F^^ : F^^^ : F^^ ~ 1 : 1 : 1 would 
be modulated to F^^ : F^^ : F^^ ~ 1 : 1 : 1 reaching Earth because the wavelength for i/^ — > i^r oscillation is much 
smaller than that of the oscillation solution associated with the atmospheric solar data. A detailed analysis of the 
detector signatures for various oscillation scenarios is presented in Refs. and ^T\ . 

Other extensions to the Standard Model could influence the predicted spectra. The existence of a supersymmetric 
(SUSY) sector would contribute extra decay channels in the fragmentation of the X-particle and jets. (See |Q and 
references therein for possible approaches to modeling SUSY channels for application to UHE decays). The net effect 
would be a skewing of the primary flux spectra to lower energies and a decrease in the energy reprocessed in high 
energy cascades. For example, in the p=l TD model studied in Ref. |Q with mx = 10^^ GeV, including a SUSY 
sector decreases the spectra of all species expected at Earth by an order of magnitude at _E > 10~^mx, shifts the 
UHE peaks in the spectra to lower energies by at least an order of magnitude and increases the amplitude of the 
Ue, Vfi, p and p spectra by about an order of magnitude at energies below the peaks. The predicted E ^ 10^^ eV 
photon spectrum, predominantly generated by EM cascades, though is unchanged because the same fraction of initial 
X-particle energy is reprocessed into low energies in the EM cascades with or without a SUSY sector in that model. 
In the case of the VHS scenario, these remarks must be convolved with the suppression of the UHE 7-ray and nucleon 
spectra due to the nearest string being farther than the relevant UHE attenuation lengths. 

To date there is no experimental evidence for a SUSY sector and no unique theoretical predictions for the architecture 
of the SUSY sector and superparticle properties such as mass. Other extensions to the Standard Model, for example 
Technicolor or Leptoquarks, which also increase the number of degrees of particle freedom would be expected to 
similarly skew the spectra to lower energies. As well, extensions to the Standard Model may modify the high-energy 
cross-sections governing the interactions of the emission with the cosmic ambient media. 




(35) 
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V. CONCLUSIONS 

We have computed the ultra-high energy 7— ray, neutrino and cosmic ray fluxes in the VHS cosmic string scenario, 
in which the long string network loses its energy directly by particle emission. This is in contrast to the standard 
cosmic string model with a scale- invariant distribution of string loops which lose energy predominantly by gravitational 
radiation. The predicted particle fluxes are much larger in the VHS scenario. 

The predictions for the particle fluxes depend on extrapolating the QCD fragmentation function to energies much 
higher than current accelerator energies. New physics between the QCD scale and the scale 77 = at which 
the defects are formed could significantly affect the final particle distributions. We parametrize this uncertainty by 
introducing a scale mj as the energy scale to which the QCD fragmentation functions can be extrapolated and regard 
TO J as the initial jet energy. The smaller the value of toj, the larger the number of jets which are generated by the 
initial emission from the cosmic string. We calculate the resulting particle fluxes as a function of both mj and /i. A 
new aspect of our work is the computation of the significant tau neutrino fluxes directly produced by decay in the 
particle jets. 

Our calculations show that the predicted flux F{E) of 7— rays, electron and muon neutrinos and cosmic rays scales 
as as in other defect models, whereas the present observations show a much weaker increase with 

E. Hence, the most stringent constraints on topological defect models come from the highest energies for which data 
exist. 

Setting Toj — the VHS strings with G/x > 10^^*^ produce electron and muon neutrinos in excess of the UHE 

observations. Thus, assuming that the strings evolve as in the VHS scenario, models with Gfi > 10^^" or, equivalently 
a symmetry breaking scale of ry > 10^'^ eV, are ruled out by the UHE data. (This result is in accordance with the limit 

on Gfi in Ref. (2^). At a given particle energy, the predicted fluxes reaching the Earth scale as rrij^^'^ and increase 
with G/i. Lowering to,/ increases the disagreement between predictions and observations, although the upper cutoff 
on the predicted spectra diminishes. 

A consistent but slightly stronger constraint on the VHS string scenario comes from the cascading of the electro- 
magnetic emission off the cosmic radiation backgrounds and the potential conflict with the observed diffuse 3 x 10^ 
eV ^ E ^ 10^^ eV gamma ray background as probed by EGRET. This requires Gfi ^ 10~^^. 

We conclude that, generically, GUT-scale strings are ruled out in the VHS scenario. However, VHS models with 
lighter Gfj,, corresponding to lower symmetry-breaking scales, may be contributing interestingly to the observed UHE 
events. We reiterate that the VHS scenario is not universally accepted as giving the correct dynamics of cosmic 
strings. 

Finally, we note that the uncertainties inherent in the calculation of the expected fluxes in any TD model mean 
that the absolute and relative flux in each particle species should not be regarded as accurate to more than an order 
of magnitude in any model. The list of uncertainties includes: the mass and other properties of the X-particle(s); 
the form of the potential of the scalar field at string formation; the branching ratios and decay channels of the X- 
particle(s); the validity of the extrapolated fragmentation functions; the possibility of new particles (e.g. a SUSY 
sector) and other extensions to the Standard Model which may influence particle production, decay and propagation; 
the strength of the extragalactic magnetic field; the intensity of the cosmic IR background; and the Hubble constant. 
Taken together, these factors have the capability to increase or decrease the absolute and relative fluxes expected at 
Earth. 
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NEUTRINO FLUX FROM MUON DECAY 

In the main text we computed the flux of particles resulting from two-body decay of charged pions in a jet. The 
muons which are produced in this decay are unstable and in turn decay via a three body decay process producing 
electrons, muon neutrinos and anti-electron neutrinos. In this Appendix we compute the full spectra neutrinos from 
the TT^ ^ fi^ ^ v decay chain. 
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Because of the finite muon rest mass, the muon and muon neutrinos created in the first stage of charged pion decay 
cannot carry an arbitrary fraction of the pion energy. In the limit E.^ 3> m^^, the energies of the decay products he 
between 

E e [0,a-^E^] (36) 

for muon neutrinos and 

E e [rE^^E^] (37) 
for muons |3^, |3^. Here, r — (^^)^ = 1 — where a is given by (^. The muon distribution is then 

= aE;^ ■ (38) 



dE^ 



In order to obtain the i^^-flux, we integrate (5|) over the range (|36[), resulting in (|^). To obtain the flux of muon 
decay products, we follow the method of Ref. [39|. In the laboratory frame the distribution of Icptons with energy 
E = yEf^ produced by the decay of a muon of energy i?^ is given by [Gaisser p. 92] 



dn 
dy 



— =9o{y)-PMy) (39) 



in the limit i?„ 3> m„,TOe, where 



for and 



ffo(y) = ^ - 3y2 + ^y^ 



1 8 

9i{y) = o -3y'+ 



giiy) = -2 + 12y - 18y^ + 8y^ 



for Ve- In (|3S|), is the projection of the muon spin in the muon rest frame along the direction of the muon velocity 
in the laboratory frame 



P, 



2E^r 1 + r 



E^{l-r) 



Convolving the distribution (^9|) with the distribution of muons from pion decay ( pSD and integrating over muon 
and pion energies, the neutrino spectrum in the laboratory frame is 

Ernin = min[rE^,E^] . 
For simplicity we will first evaluate ( |40| ) for an initial pion distribution of the form 

dN 



Changing the order of integration. 



dE. - ^^'^^ ' « > « 



dE^ / dE^^ dE^, / dE^ , 
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the result is 



where 



dEy " a(l-r) 



1 -r 



1 + r - 



2ar /I - r" 



a - 1 V 1 - r"^ 



(41) 



/o 



2(a + 5) 5 fE^y 3 / -B, 



a+2 



a{a + 2){a + 3) 3a\nij J ' {a + 2)\mjj 3{a + 3)\mj 



a+3 



for and 



/i = 



2(1 - a) 



a{a + 2){a + 'S) 



3a \mj 



3 fE, 
{a + 2)\mj 



a+2 



3(a + 3) \ mj 



a+3 



12 



2 



a(a + 2)(a + 3) a \mj 



6 fE, 
(a + 2) l^mj 



a+2 



E, 



(a + 3) \r7ij 



a+3 



12(a- 1) 



2fE, 



a(a + l)(a + 2)(a + 3) a\mj 

a+2 , 



+ 



18 



E,, 



(a + 2) \mj 



8 fE, 
(a + 3) \mj 



12 

(a + 1) Vmj 

a+3 



(42) 



for v^. Because muons and antimuons are created with opposite spin polarization in charged pion decays, this spectrum 
also applies to the antiparticle decay chain. 

The distribution of pions produced in QCD jet decay can be approximated by the polynomial fragmentation function 
(^. Evaluating (^l|) for the appropriate values of K and a and extending the analysis to a < 0, we find the final 
distribution of neutrinos produced by jet decay to be approximately 



dN, 



dE, 



8mj 
- 7.16 



12.48 + 0.44 



mj 



6.12 



-1/2 



Eu 

mj 



1/2 



E^.. < 



mj 



where a is given in (Q) , for the muon neutrinos produced in the first stage of the pion decay, 



dE^.. 



- 7.33 



11.83 + 0.48 

E, 



E^ 
mj 



-3/2 / ^ 

-5.80 



-1/2 



0.15 



mj 

Eu 
mj 



1/2 ■ 

+ 0.97 

\m,j 
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Ei,^^ < mj 



for the muon neutrinos produced by the decay of muons in the pion decay, and 



dE^^ 



8mj 



12.68 + 0.47 



.30 



1.67 



E^ 
mj 

Eu 
mj 



1/2 



E^ 
mj 



0.11 



- 0.45 



-3/2 

E^ 
mj 

E, 
mj 



5.96 



mj 



-1/2 



Ev^ < mj 
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for the final electron neutrinos. 
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